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Isospin equilibration in multi-fragmentation processes is studied for the system '^^Cl +'^^ Si at 

00 ^ 

. 40 MeV/nucleon. The investigation is performed through semiclassical microscopic many-body 



calculations based on the CoMD-II model. The study has been developed to describe isospin equi- 
libration processes involving the gas and liquid "phases" of the total system formed in the collision 
processes. The investigation of the behavior of this observable in terms of the repulsive/attractive 
action of the symmetry term, highlights many-body correlations which are absent in semiclassical 
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An interesting subject related to the Heavy Ions Isospin physics [ll is the process leading 
to the equilibration of the charge/mass ratio between the main partners of the reaction. The 
so called "isospin diffusion" phenomenon has been indicated as the associated mechanism 
acting between the reaction partners 

mm 

. In particular, in basically binary processes 
occurring after the collision of Sn 124 and 112 isotopes [3(], evidence of partial equilibrium in 
the charge/mass ratios of the quasi-projectile (QP) and quasi-target (QT) has been deduced 
through the study of the iso-scaling parameters related to the isotopic distributions. In this 
case dynamical calculations based on the BUU model jo] show that the degree of equilibration 
depends on the behavior of the symmetry potential as a function of the density. The 
analysis in this kind of studies, however, is based on the linear relation between the iso-scaling 
parameter and the relative neutron excess Y of the emitting sources (an assumption typical of 
several statistical models), moreover, it is also assumed that both quantities weakly depend 
on the secondary statistical decay processes. However, in general, this last condition depends 
on the fragment excitation energies, or temperatures, and on the distinctive features of the 
models used to simulate the second stage decay [7]. In this work we want to describe the 
isospin equilibration process, looking at the whole system, by using the dynamical variable 
'^{t) = J2i=i 'v'i- The sum on the index i is performed on all the protons of the system. 
17 (t) corresponds to the time derivative of the total dipole of the system. The velocities IjI 
are computed in the center of mass (cm.) reference system. Several studies were based on 
this dynamical variable to describe pre-equilibrium 7-ray emission (see Refs.flSflQ and 



references therein ). Various reasons dictate this choice to describe also isospin equilibration 
processes. 

- i) After the pre-equilibrium stage at the time tp^e, when a second stage characterized by 
an average isotropic emission of the secondary sources (statistical equilibrium) takes place, 
the ensemble average satisfies the following relation: 

^{tpre) = V(t > Ve) = 1^ . The 

average value of the dynamical variable at tpre is in fact invariant with respect to statistical 
processes and therefore is an interesting observable to be investigated. In particular, 
l7 can be expressed as function of the charge Z, mass A, average multiplicity rfiz^A and 
the mean momentum {P)z,A of the detected particles having charge Z and mass A in the 
generic event: 

Z 

Z,A 
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^z,A ^ rnz,A{?)z,A ^2) 

'W)z,A^Z,A 

is the correlation function between the multiplicity and the mean momentum. This 

{P) 

correlation function plays a determinant role for the invariance property and therefore asks 
for an analysis event by event in which many-body correlations can not be neglected. 

W , for symmetry reasons, he on the reaction plane. As suggested from Eq.(l), it is 
directly linked with a weighted mean of the charge/mass ratio. It takes into account also 

the average isospin flow direction through the momenta {P) z A' 



-ii) In the general case, we find attractive the following decomposition: = V^g -|- 
where with \^ a and V^l we indicate the dipolar signals associated to the gas 
"phase" (light charged particles) and to the "hquid" part, corresponding to the motion of the 
ensemble of the produced fragments (QP,QT, if any, and intermediate mass fragments). The 
signal V^GL is instead associated to the relative motion of the two "phases". By supposing, 
for simplicity, that the gas "phase" is formed by neutrons and protons, l7 can be further 
decomposed as: 



V = ^ V ^ H ^ V cm,LG + Vr,L (3) 

In the above expression the first term represents the contribution related to the neutron- 
proton relative motion of the gas "phase" expressed through the relative velocity if^^, 
the second term is related to the relative motion lfcm,LG between the centers of mass of 
the "liquid" complex and the "gas"; the last term represents the contribution produced by 
the relative motion of the fragments. A similar expression can be obtained including in 
the gas "phase" other light particles. From this decomposition we can see how the isospin 
equilibration condition = 0), for the total system, requires a very delicate balance 
which depends on the average neutron excess of the produced "liquid drops" Yl, on the one 
associated to the gas "phase" Yq and on the relative velocities between the different parts. To 
enlighten the role played by some of the terms reported in Eq.(3), we can discuss the idealized 
decay of a charge/mass asymmetric source through neutrons and protons emission (or the 
case in which the liquid drops are produced through a statistical mechanism \^r,L — 0). 
Moreover, for simplicity, we can consider uncorrelated fluctuations between the velocities 
and neutron excesses. In absence of pre-equilibrium emission or for identical colliding nuclei 



( V cm,LG — 0) the second term of Eq.(3) is zero, and the isospin equilibration will require a 



proton-neutron symmetric gas "phase" and/or absence of relative neutron-proton motion. 



For non identical colliding nuclei, if pre-equilibrium emission is present, then v cm,LG 7^ 0. In 
this case, if Yq 7^ Yl, as due to the isospin "distillation" phenomenon, the first term has to be 
necessarily different from zero and it will contribute to the neutron-proton differential flow 
Fnp (see the following). Therefore, for the whole system, the isospin equilibration process 
can not be explained only in terms of isospin diffusion and drift between the main partners 
of the reaction. It is in fact necessary to take into account the gas "phase" contribution 
which we can regard as a kind of "dissipation" with respect to the system formed by the 
liquid part. In particular in this work, as an example, we will discuss the results obtained 
through the Constrained Molecular Dynamics-II approach (CoMD-II) 1^, 13| applied to the 



charge/mass asymmetric system -|-^^ Si at 40 MeV/nucleon. The study is performed 
by using different options for the symmetry potential term U"^ . Nevertheless, to clearly 
understand the dynamics of the investigated system, it is necessary to dedicate the following 
part to a discussion of this term. In particular, we will illustrate that the existence of many- 
body correlations, as produced in our approach, strongly affects this part of the interaction. 

In the framework of the present version of the CoMD-II model the isospin momentum 



dependent part of the total energy is included through the Pau 
the potential symmetry term has been implemented as follows 



i principle constraint and 



= ^F'{s)Pm (4) 

« ^NN I ^PP ^NP ^PN /r\ 

PM = P+P-P-P [0) 

iCKJCK 

p""""' = E K,K' = N or P (6) 

F\s) = SUffl (7) 

F'{s) = 1 SUff2 (8) 

F'is) = (— Soft (9) 

Sg.s. 

^ iJCKUK' 

^ = TTT E P.. (10) 

With the superscript N and P we indicate the neutron and the proton ensemble of nucleons 
respectively and A is the total mass of the system, pij represents the normalized Gaussian 

overlap matrix elements 12] typical of the quantum molecular dynamics approach The 



label g.s. indicates the ground state configuration. The positive form factors F'{s) are 



such that F'{s)j^ has the same functional form hke F{u) proposed in 15| (see Eq.(3) 
in Ref ■ [isl ) . They have been used extensively in mean-field (M.F.) approaches. Apart 
from a normalization factor, equal to |, s is related to the average (with respect to the 
number of two-body interactions) overlap integral . For the system under study and for 
compact configurations s is well approximated by the average one-body density p. The 
factor (3m instead arises naturally from the many-body approach. It takes into account 
explicitly, through the last two negative terms in Eq.(5), that the microscopic two-body 
nucleon-nucleon interaction in the isospin singlet states is more attractive than the one 

n 

related to triplet states p] (ctsym = 72MeV). For the moderately asymmetric system, here 
investigated, self-consistent calculations, including the ones for the searching of the g.s. 
properties, produce a negative value of Pm- This result is due to the combined actions 
of the more attractive force, for the isospin singlet states, and the repulsive action of the 
Coulomb interaction for proton-proton couples. We have verified also that the Pauli principle 
plays a role even if with a smaller extent. These correlations, in fact, tend to increase the 
average (with respect the number of two-body interactions) neutron-proton overlap integral 
p^^ and to decrease the overlap p^^ and p^^ related to the neutron-neutron and proton- 
proton couples. These average overlap integrals are defined as: p^^' = p^^' /II' with / and 
/' equal to and/or P. In our calculations we have verified that p'^^ = = a— — ^ — - 
with a ^ 1.1 — 1.2. Therefore, at low asymmetry, it is enough a small correlation effect 
to produce a negative value of Pm, as due to the structure of Eq.(5). However, it results 
that Pm decreases by increasing the symmetry of the system (the number of singlet states 
increases) and, on the contrary, it increases by decreasing the symmetry of the system. These 
changes follow a parabolic dependence with respect the x = ^ ~ ^ variable. This means 
that Pm contains, independently from the sign, the right behavior necessary to explain the 
differences in the binding energies of isobars nuclei and it is able to generate the so called 
"isospin distillation" phenomenon js] in the dynamics of hot sources. However, the negative 
sign of Pm has relevant consequences in to determine the effects of the different options 
for U"^ on the dynamics of the investigated processes. For example, contrary to the semi- 
classical M.F. description, the so called Stiff cases (in spite of this change of sign, we retain 
the same nomenclature to indicate the different options) show a decreasing behavior with the 
density. The Soft case, here investigated, even if negative, maintains the increasing pattern 
typical of a M.F. approach (see Fig.5(c)). We conclude this part by briefly describing how in 
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M.F. and/or in Liquid Drop Model these correlations are washed out. For this purpose we 
assume, for simplicity, A,N, Z 1, and we rewrite Pm as: Pm = N'^p'^^ +Z^p^^—2NZp^^ . 
M.F. approximation is now easily obtained by supposing p^^ = — '^^^2^ ''^ = P = ^P be 
independent on I and I'. This means that the correlation producing the aforementioned 
differences in the matrix elements p^^' are averaged through the overlap integral p, so that 
■ Pm Ja(-^~^)'^ ~ ^"'i^ , where p„, pp and p represent the neutron, proton 

and total densities respectively. V is the volume. In this case, therefore, we obtain that the 
interaction depends locally only through one-body densities and it is positive defined. Taking 
into account that for the present system, in compact configurations, s = p, by substituting 
the right end of the above relation in Eq.(4) and dividing by the volume V, we get a 
precise correspondence with the symmetry energy density used in Ref. [isl if Cq = 27MeV. 
This value includes the reduction of the Fermi motion as due to finite size effects and 
corresponds to an Sq value (see Eq.(4) of Ref. IS]) of about 36 MeV. Finally we note that 
these correlations not only change the sign of the symmetry energy, but they heavily affect 
also the strength. In particular, according to Eqs.(4-10) and to the definition of Pm in 
terms of the average overlap integral p^^' it is straightforward to evaluate the contribution 
of Pm by assuming 10% of correlation, i.e. a = 1.1. By considering a total system like 
the one investigated in the present work (A=68 Y=0.088) at normal density we obtain 
Pm — — 0.43/m~^ in the correlated case and Pm = 0.065fm~^ in the M.F case. According 
to this evaluation we observe that also with a correlation value of the order of 1% the 
induced effect on the symmetry term is not negligible. This strongly suggests that the 
qualitative effects discussed in the present work can be still present by using other kinds of 
effective interactions. We have to note also that, according to the increasing behavior of the 
symmetry energy with the charge/mass asymmetry, for systems with mass around 64 units, 
the Pm factor becomes positive for relevant asymmetry \Yc\ > 0.25 restoring the repulsive 
behavior of the asy-stiff cases here investigated. 

Now we can discuss the results concerning the isospin equilibration process for the system 
'^^Cl +'^^ Si at 40 MeV/nucleon. For the system under study, in Fig. 1(a) and Fig. 1(b) 
we show the average total dipolar signals as evaluated through CoMD-II calculations along 
the z beam direction V'^ and along the impact parameter direction x, , respectively. The 
reference frame is the cm. one. The impact parameter b is equal to 4 fm. 
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FIG. 1: Average dipolar signals along the z and x directions are plotted as a function of time in 
the intervals 0-150 fm/c (panels (a) and (b)) and 100-400 fm/c. The three kinds of lines refer to 
the different options (see eqs.(7-8)) describing the symmetry potential. 

In these figures the average dipolar signals are shown for the first 150 fm/c. Different 
lines refer to different symmetry potentials. The isospin independent compressibility is equal 



to 220 MeV according to Ref. 12] ■ In the first 150 fm/c we can see that in all the cases 



wide oscillations are present. They are responsible for the pre-equilibrium 7-rays emission 







9|. The damped oscillations tend to smaller and constant values (within the uncertainty 
of the statistics of the ensemble average procedure) as can be seen in Fig. 1(c) and Fig. 
1(d) in which the dynamical evolution is followed from 100 fm/c up to 400 fm/c. The time 
interval in which the stationary behavior is reached corresponds to the average time for the 
formation of the main fragments. In Fig. 2(a), for the different symmetry potentials, we plot 
the charge Z distributions. The Stiffl and Stiff2 options clearly show a multi-fragmentation 
pattern, while for the Soft case the behavior is closer to an exponential trend making reliable 
a "vaporization" scenario. 

We now briefly comment the results already shown in Figs. 1(c) and 1(d). We can see 
that the asymptotic values sensitively depend on the options used for the symmetry term. 
For example, concerning the dipolar signal along the z direction, the difference between the 
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FIG. 2: In the panel (a) we show the charge Z distributions of the reaction products for the 
different options describing the symmetry potential (see eqs.(7-8)). In panels (b) and (c), for the 
Stiffl and Soft cases, the z components of the dipolar signals are shown for the liquid and gas 
"phases". 

Stiffl and the Soft cases is about 85% while between the StifFl and the Stiff2 options is 
about 53%. The changes of the dipolar signals along the x direction are smaller of the order 
of 10-30%. To understand the role of the "gas" particles, in Fig. 2(b) and in Fig. 2(c) 
we show with open circles and with closed circles the values of dipolar signal along the z 
direction, V^q (gas "phase") as a function of time and V\ (hquid "phase"), for the Stiffl and 
Soft parameterizations respectively. We can see that the contributions related to the "gas" 
are negative (light particles coming essentially from the target) . They are quite effective in 
to reduce the absolute value of the total dipolar signals driving the total system through the 
isospin equilibration. It is also possible to sec that shows a larger absolute value for the 
Soft case with respect the Stiffl option. This is due to the more repulsive effect obtained in 
the Soft case with respect to the other ones. In Fig. 3(a) we show, as a function of time, 
the average overlap integral s for the three cases. 

In Fig. 3(b) and Fig. 3(c) we show the corresponding total average potential C/*°* and 
the symmetry one U^. From Fig. 3(a) we can see that the maximum value of s, Smaxj is 
reached in about 30 fm/c. The Soft case shows the lower s^ax value, the Stiffl the higher one 
(about 1.25 times the Sg,a. value). Accordingly, f/*°* displays the less attractive behavior for 
the Soft case (see Fig. 3(b)). The reasons of these differences can be understood by looking 
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FIG. 3: (a) average overlap integral s, (b) average total potential, (c) average symmetry potential 
are shown as a function of time for the different options describing the symmetry potential (see 
eqs.(7-8)). 



at Fig. 3(c). f/^ shows, in fact, a clear repulsive behavior as a function of the density for 
the Soft case while the Stiff 1 case gives rise to an attractive effect. At later times, when 
the density is decreasing and the fragment formation process takes place for s < Sg,s,, the 
behavior of 17^ in the different cases is reversed. The Stiff2 case shows always an attractive 
behavior with intermediate values concerning the strength. The differences observed around 
30 fm/c play a crucial role in to determine the later evolution of the system 161]. This can 
be understood by looking at the figures 2(a) and 2(b). In these figures it is in fact clearly 
visible that the dipolar signals associated to the gas "phase" start to develop just at 30 
fm/c. The less attractive behavior f/*°* for the Soft case determines a lower collision rate 
and a higher multiplicity of light particles. For example, the total number of neutrons and 
protons emitted at 160 fm/c is, on average, about 15 for the Soft case and 11 for the Stiffl 
one. The gas "phase" is more rich in neutron for the Soft case {Yq = 0.23) with respect to 
the Stiffl option {Yq = 0.18). These values are rather large as compared with the neutron 
excess fraction of the total system and with the one related to the biggest fragment. This 
clearly indicates the occurrence of the so called "isospin distillation". As above mentioned 
the dipolar signals associated to the relative neutron-proton motion for the gas "phase" can 
affect the neutron-proton differential flow Fnp As an example in Fig. 4(a) for the Soft 
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FIG. 4: (a) Differential neutron-proton flow F„„ as a function of (— ^ — )cm (see text), (b) The 

^ Vbeam 

same quantity is plotted after the corrections given by the relative neutron-proton motion related 
to the gas "phase", (c) For different symmetry potentials we show the relative changes of the 
average total dipolar signal along the z component as a function of the impact parameter. 

and Stiffl cases we show the neutron-proton differential flow as a function of the particles 
rapidity y normalized to the projectile one Ubeam- 

The rapidity values arc evaluated in the cm. reference system. In Fig. 4(b), we show 
the same quantities corrected for the average free neutron-proton relative motion. We can 
see that the differences are considerable and, in this case, they are due essentially to the 
X component of the relative motion. We conclude this study by showing in Fig. 4(c), for 
the different symmetry terms, the relative changes of the average total dipolar signals 
along the z component as function of the impact parameter. It is clearly visible the higher 
sensitivity obtained for the impact parameter range h = 3.5 — 5.5 fm. CoMD-II calculations 
show that, for more central collisions, even if the relative changes in the overlap integrals 
are slightly more pronounced, the increasing of the collision rate produces a more damped 
mechanism with respect to the dipolar degree of freedom. This determines a smaller sensi- 
tivity to the change of the symmetry energy. On the contrary, for more peripheral reactions, 
the dipolar signal shows a smaller damping and a related smaller degree of equilibration. In 
this case, however, the changes in the relative overlap integral s are smaller so that the as- 
sociated changes in U'^ are less pronounced. These suppression mechanisms and the related 
existence of an optimal range of impact parameters for an enhanced visibility of the isospin 
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effects seem to us quite general, and clearly depending also on the way in which the collision 



term is treated 



13|. 



In summary, in this work the isospin equilibration process for the asymmetric charge/mass 
system "^^Cl Si at 40 MeV/nucleon has been investigated by studying the ensemble aver- 
age of the time derivative of the total dipole \^ as evaluated through CoMD-II calculations. 
Some general properties of this quantity have been discussed. In particular, it allows to give 
a consistent definition of isospin equilibration also in complex reactions evolving through 
multi-fragmentation processes. CoMD-II calculations show that the asymptotic values of 

for these processes are quite sensitive to different symmetry potential options; moreover, 
for this system, the dipolar contribution associated to the pre-equilibrium charged particles 
emission is relevant in to determine the value of V^. Detailed CoMD-II calculations also 
enlighten fundamental differences in the behavior of the symmetry potential with respect to 
a description based on a M.F. approach. The negative sign of the leading Pm factor in the 
symmetry potential represents a typical example in which the ensemble average, performed 
on many realizations of a semiclassical many-body dynamics approach, gives arise to deep 
differences with respect to the average behavior evaluated with a M.F. approach. These 
deviations are due to the highly correlated neutron-proton motion generated by the symme- 
try term and they are clearly beyond the one-body description. Moreover, the calculations 
presented in the present work suggest the existence of relevant critical charge/mass asym- 
metries starting from which the symmetry potential changes sign producing a transition 
from a density attractive to a density repulsive behavior for the asy-stiff cases. We conclude 
by observing that this subject needs further detailed studies in the next future together 
with a specific investigation about the effects of these correlations on the behavior of others 
isospin-dependent observables. 
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